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Introduction 

This  paper  is  concerned  with  the  problem  of  analytically  modeling  the 
acoustic  fields  that  control  the  directional  response  of  transducer  array 
elements  mounted  on  the  flat  face  of  an  elllpsoldally  capped  cylinder,  such 
as  that  depicted  schematically  in  Figure  1.  The  basic  objective  Is  to 
Investigate  the  extent  to  which  the  back-to-front  ratios  of  Individual 
transducer  array  element  directivity  patterns  can  be  minimized  by  reducing 
the  diffracted  fields  which  reach  them  from  the  rear  of  the  body  via 
propagation  over  the  curved  surface. 

The  measurement  procedure  to  be  modeled  is  plane  wave  insonlflcatlon  of 
the  array  elements  with  the  source  In  the  far-fleld  of  the  body,  as  depicted 
schematically  in  Figures  la  and  lb;  the  former  with  the  source  in  the  'Insonifled 
zone  of  the  array*  (l.e.,  0 £ 0 it/2)  , depicting  the  combined  incident, 

reflected,  and  diffracted  fields  normalized  with  respect  to  the  incident 
field;  the  latter,  with  the  source  In  the  'shadow  zone  of  the  array*  (l.e., 

Tr/2  < d ^ it),  depicting  the  field  diffracted  around  the  curved  surface  of 
the  body  to  the  array  element,  again  normalized  with  respect  to  i ^ incident 
field.  Note  that  when  the  terms  'insonifled  zone'  or  'shadow  zone'  of  the 
array  are  used  here,  the  field  at  the  source  is  regarded  by  reciprocity  to 
be  produced  via  radiation  from  the  array. 

In  order  to  describe  the  diffracted  field  analytically,  we  make  use  of 
1 2 

Keller's  ’ Geometrical  Theory  of  Diffraction  (GTD).  This  method  effectively 
permits  ray  tracing  over  curved  impedance  surfaces  where  the  only  part  of  a 

3 

ray  tube  that  can  propagate  are  tangent  rays  (l.e.,  Franz  waves,  analogous 
to  Stoneley  waves  which  occur  at  almost  grazing  incidence  on  insonifled  flat 
elastic  surfaces  immersed  in  liquids) . Since  the  diffracted  field  sheds 
rays  at  each  point  of  tangency  as  it  propagates  along  a curved  surface,  the 


UNCLASSIFIED 


UNCLASSIFIED 


-10-» 


November  17,  1977 
FHF;cjp 


energy  flux  diminishes  in  a manner  phenomenologically  attributable  to  an 
exponential  field  decay  with  geodesic  distance  from  the  point  of  Insonlflcatlon. 
Now  the  'attenuation  coefficient'  of  this  decay  Is  dependent  on  the  local 
radlus-of-curvature  of  the  body  and  on  the  local  surface  Impedance,  the 
latter  dependence  being  a consequence  of  the  manner  In  which  the  Impedance 
affects  re-radlatlon  of  the  diffracted  field  at  each  point  of  the  surface. 

Thus,  for  example,  the  diffracted  field  on  an  acoustically  rigid  curved 
surface  Is  re~radlated  at  each  point  of  the  surface  by  a monopole  source 
distribution.  This,  In  turn.  Is  more  efficient  than  the  dipole  source 
distribution  associated  with  an  acoustically  soft  (l.e. , pressure  release) 
surface.  Consequently,  the  'attenuation'  per  wavelength  of  surface  rays 
propagating  on  an  acoustically  soft  curved  surface  Is  considerably  greater 
than  that  on  an  acoustically  hard  curved  surface.  Since  maf.erj&ls  which 
are  'rigid'  in  air  appear  acoustically  more  compliant  in  water,  we  note  at 
this  point  that  'rigid'  In  water  Implies  an  Impedance  of  at  least  an  order- 
of-magnltude  greater  than  the  characteristic  Impedance  of  the  fluid. 

The  'diffraction'  and  'attenuation'  coefficients  of  Keller's^ 

Theory  are  obtained  from  'canonical  problems'  l.e.,  those  that  are  amenable 

to  exact  solution,  such  as  creeping  wave  solutions  for  separable  coordinate 

surfaces  of  the  Helmholtz  equation.  In  particular  surfaces  of  constant 

curvature  ( e.g.  Infinite  cylinders  and  spheres).  For  these  basic  or 

canonical  surfaces  the  creeping  wave  solutions  are  obtained  by  re-expresslng 

the  exact  eigenfunction  solutions  (whose  rate  of  convergence  deteriorates  with 

Increasing  values  of  kp2  In  the  form  of  Luneberg-Kllne  series  whose  term^ 

decrease  as  (1/kp)”,  n » 1,2,  . . . Whenever  kp.»  1 , It  follows  that 
* * 

only  the  first  few  terms  of  these  reconditioned  series  are  required  to 
provide  an  accurate  representation  of  the  diffracted  field.  In  the  case 
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canonical  series  can  also  be  used  to  describe  the  diffracted  field  over 


local  regions  of  the  surface  where  the  radlus-of-curvature  Is  slowly  varying 


although  a note  of  caution  should  be  added  In  the  light  of  Lepplngton's 


analysis.  Generally  speaking,  however 


In  terms  of  Its  Impedance  and  differential  geometry  as  shown  by  Levy  and 


Before  concluding  these  remarks  a point  of  confusion  between  electro' 


magnetic  and  acoustic  analogues  should  be  mentioned.  This  point  concerns 


the  application  of  the  electromagnetic  Impedance  boundary  conditions  to 


acoustic  problems,  l.e.  as  given  by  Keller 


where  (|i  Is  the  scalar  e.m.  potential,  and  Z Is  the  e.m.  surface 


Impedance  relative  to  the  characteristic  Impedance  of  the  surrounding 


medium.  In  an  acoustic  field,  however,  since  the  particle  velocity  Is 


given  by  definition  as  v > 7(^  and  for  sinusoidal  disturbances  the 


pressure  Is  given  by  p ~ -Ip  c k(^  , the  Impedance  relationship  p « z |v 


where  (P  Is  the  acoustic  field  potential,  and  Z * z /pc  Is  the 


acoustic  surface  Impedance  relative  to  the  characteristic  Impedance  of 


the  fluid  medium.  From  Inspection  of  Equations  (la)  and  (lb)  It  follows 
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that  whenever  the  acoustic  analogue  of  a corresponding  e.m.  diffraction 
problem  Is  Invoked,  Z should  be  replaced  by  Z ^ yielding  the  relation- 
ships  between  diffraction  and  attenuation  coefficients  summarized  by  Pathak 
and  Kouyoumj Ian. ^ 

Theory 

We  begin  by  considering  the  axlsymmetrlc  Insonlf Icatlon  of  an 
elllpsoldally  capped  cylinder  such  as  that  depicted  schematically  In  Figures 
la  and  lb  by  a point  source  In  the  far-fleld  of  the  body.  For  the  case  of 
a circular  piston  of  diameter  d mounted  on  the  flat  front-end  surface, 
the  received  field  can  be  expressed  as, 

P " ?!  + Pr  + Pd  » C2) 

where  Pj^>P^>  P^  denote  the  Incident,  reflected,  and  diffracted 

fields,  respectively.  When  the  source  Is  located  In  the  first  quadrant 
(l.e..  In  the  Insonlfled  zone  of  the  array,  0 0 £ — ) , If  the  length  of 

the  body  Is  sufficient  to  ensure  the  effective  attenuation  of  diffracted 
waves  before  they  arrive  at  the  receiver  elements  via  propagation  around 
the  body,  the  field  received  by  an  element  from  0°  to  90”  will  only 
depend  upon  the  Incident  and  reflected  waves.  In  this  Instance  therefore, 

P “ P^  + Py  f'  Pd  " ® ’ 0 1 ® 1 f • 

On  the  other  hand,  when  the  source  Is  In  the  second  quadrant  (l.e., 

^ < 6 £ tt) , only  diffracted  waves  reach  the  array  elements,  so  that  In  the 
shadow  zone  of  the  array. 
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P - 


TT 


j < 0 £ ir 


(3b) 


In  the  third  quadrant  the  field  is  again  given  by  Equation  (3b)  and  in  the 
fourth  quadrant  by  Equation  (3a) . 

Between  the  Insonlfled  zone  and  the  deep  shadow  zone  of  the  array  there 
is  a region  described  as  the  transition  zone  which  approximately  spans  the 
sector  "kH  - (kp^°/2)  £ 0 £ '*f/2  + (kp^°/2)  . In  this  region,  the 

field  is  governed  by  Incident,  reflected,  and  diffracted  waves.  One  con- 
sequence of  the  transition  zone  is  that  the  directivity  pattern  or  directional 
response  of  a receiver  element  normalized  with  respect  to  its  response  on 


boresight,  is  reduced  by  50Z  at  9 - ^72  on  uniformly  hard  (l.e.,  Z > <») 

a 


and  soft  baffles  (l.e.,  Z >0).  Under  these  conditions,  as  Pathak  and 

a 


Kouyoum j Ian”  have  shown  for  analogous  e.m.  problems,  the  transition  zone 
field  can  be  described  analytically  by  means  of  Fock  functions.  However, 
the  derivation  of  a suitable  transition  zone  approximation  for  an  arbitrary 
impedance  baffle  has  yet  to  be  resolved,  and  consequently  the  exact  value 
of  a receiver  element ' s response  at  6 - it/2  remains  unspeclfable.  For 


this  reason,  in  our  subsequent  numerical  analysis,  we  compute  the  response 


from  0 to  Tr/2  - (kp^°/2)"^^^  and  from  Tr/2  + (kp^°/2)"^'^^  to  n , 


Joining  the  two  regions  graphically  via  matching  slopes.  Although  this 
procedure  cannot  be  rigorously  justified,  it  does  not  in  any  way  detract 
from  our  computation  of  the  response  in  the  deep  shadow  zone,  which  is  the 
primary  task  of  this  investigation. 

Concentrating  on  the  field  in  the  insonlfled  zone.  Equation  (3c)  can 
be  re-expressed  in  the  form. 


P ■ Pj^d  + R(0)}  . 


0 < 0 < T ^ 


(4) 
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where  the  reflection  coefficient  R(6)  is  defined  as. 


Z (0)  cos (6  “ 0 ) - 1 

“ z’’(0)  cos(0  - 0°)  + 1 ’ 
P o 


Z (0)  being  the  angle-dependent  plate  impedance  of  the  flat  face, 

P 

normalized  with  respect  to  the  characteristic  impedance  pc  of  the 
medium  in  which  the  body  is  located,  i.e.,  as  summarized  by  Hayek  and 
Stuart,^ 


Z^(0)  - -l(k  h)(p  c /pc  ){n^  ,in^(0  - e ) - 1}  , k^  - u/c  . (6) 

P P P P o P P 


In  this  notation,  p is  the  mass  density  of  the  plate  and  c is  the 

P P 

phase  velocity  of  plate  waves  which  in  turn  is  a function  of  the  Young's 
modulus  and  Poisson's  ratio  of  the  plate  given  by 


(1  - a ")p 

P P 


Again  n ■ u/o) 


(7) 


(8) 


where 


h c 


0)  ■ — : — is  the  coincidence  frequency  for  a plate  (9) 


p of  thickness  h 


Substituting  Equation  (5)  in  Equation  (4)  and  multiplying  by  the  directional 

response  function  of  a circular  piston  of  diameter  d , the  resulting 

directivity  pattern  of  a receiver  located  at  the  point  (r  ,0  ) on  the  flat 

o o 

face  of  the  body  can  be  expressed  as. 
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D(0)  - p/p. 


2Z  (0)  co8(0  - 0 ) 


Z (0)  cos (6  - 9 
P o 


^ ^ 8in(0  - 0 ) J 

L O 


ikr  cos(0  - 0 ) 
o o 


0 < 0 - 0^  < f (10) 

— 0 — 2 


We  now  consider  the  back  response  of  a receiver  element  due  to  diffracted 
waves  which  occur  when  the  source  Is  located  In  the  shadow  zone  of  the 
array,  ir/Z  < 0 j<  tt  . From  Levy  and  Keller's^  fundamental  analysis,  this 
field  can  be  represented  via  the  Geometrical  Theory  of  Diffraction  as, 

o . I . jjLK  lk(t  + S)-[p2 

Pd  " liT  tsa^Cp^  + s)  / Jq  °n  ^^1^  °n  ^^2^®  ^^1  “ 

(11) 

The  parameters  a , a.,  p, , S,  t,  and  R which  appear  In  this  equation 
oil 

are  shown  In  Figure  2a  for  diffracted  waves  traveling  over  an  arbitrary 
curved  surface  from  the  point  of  Insonlf Icatlon  to  the  point  of  departure 

^2  • Analytical  expressions  for  the  diffraction  and  attenuation  coefficients 

Tjr  LK  5 

D and  a , respectively,  deduced  by  Levy  and  Keller  are  compared  In 
n n 

Table  la  with  similar  expressions  resulting  from  a different  definition  of 
the  Airy  Function  as  summarized  by  Pathak  and  Kouyoumjlan^  - numerical 
comparisons  being  presented  In  Table  lb.  A brief  synopsis  of  the  steps 
Inherent  In  the  derivation  of  these  coefficients  Is  also  Included  for  the 
benefit  of  the  reader  In  Appendix  A. 

For  the  case  of  the  axlsymmetrlc  elllpsoldally  capped  cylinder  depicted 
In  Figure  2b,  the  variable  radlus-of-curvature  p^  of  the  elliptical 
section  Is  given  by. 
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Flg«  2a  GTD  Parameters  of  a Diffracted 

Wave  Field  on  an  Arbitrary  Surface 
of  Revolution 


Fig.  2b  GTD  Parameters  of  a Diffracted  Wave  Field  on  an  Elllpsoldally  Capped  Cylinder 
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Equation  (B2)  ; 
n 

1/3 

3 ' q^  , from  Equation  (A20) ; 


from  Equatlon(A18) 
from  Equatlon(A19) 


Parameters 

Hard  Acoustic  Surface,  " 

Soft  Acoustic  Surface,  Z - 0 

a 

To 

0.8086166 

1.8557571 

Tl 

2.5780962 

3.2446076 

LK 

**0 

1.46935 

3.37213 

LK 

’l 

4.68417 

5.89584 

*10 

1.01879 

2.33811 

3.24820 

4.08795 

A / LK. 

1.16680 

0.0 

-0.91273 

0.0 

0.0 

-1.05905 

0.0 

1.21295 

0.53566 

0.0 

A^(-qi) 

-0.41902 

0.0 

0.0 

0.70121 

Ai(-qi) 

0.0 

-0.80311 

Table  lb. 


Numerical  Comparison  of  Parameters  Implicit  In  the  Two 
Notatlonal  Forms  of  Table  la 


I 
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+ y^3/2 


r 2 ,^3/2 

■ V } 


(12) 


O 

Following  Sachs,  who  considered  the  case  of  a source  located  on  the 
jcyllndrlcal  surface  of  the  body  rather  than  as  herein  envisaged  In  the 
far  field,  we  let 


X ■ a sin  , y ■ hj^  + b cos  . (13) 

2 

K - 1 - 2 • 

a 

Likewise,  for  far-fleld  Insonlflcatlon  of  the  body,  the  distance  t on 
the  elllpsoldally  curved  surface  between  an  arbitrary  point  of  Incidence 
and  the  exit  point  F2  of  the  diffracted  field  at  the  flat  face  Is 
given  by. 


Equation  (12)  thus  becomes. 


t 


{1  - sln^  dip  . 


(15) 


From  Table  1 and  Equation  (15),  the  Integral  of  the  attenuation  coefficients 

Tjr 

O In  Equation  (11)  over  the  Interval  t from  P,  to  P of  the 

Q X Z 

ellipsoidal  surface  can  be  expressed  as. 
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f ^ a“^(t')  dt’  - f df 

Jp,  “ Jo  “ 


■ r 


° /l  - 8in^  \|^ 


where  are  the  roots  of  the  transcendental  equation 

this  equation  being  deduced  In  Appendix  A from  the  boundary  conditions  for 
a finite  Impedance  surface.  In  practice,  however,  for  a finite  Impedance 
surface  It  Is  preferable  to  compute  the  roots  via  the  alternative  'tangent 
equation*  derived  by  Keller,^  as  outlined  In  Appendix  B. 

LK 

Now  the  diffraction  coefficients  at  the  point  of  Incidence  D (F,) 

n 1 

are  given  by. 


2'^^^  6~^^^(ab/k)^^^  sln^  (18) 

-.,2.  LK.  . LK  ,2  , LK. 

“iK<\  > ■>»  > 


I ir 

%fhllst  those  at  the  exit  point,  D (P_)  are  given  by, 

n ^ 


f^n  ,-1/2 

,,,2,  LK,  ^ LK  .2  , LK. 
%t<S,  ) + S.  *LK<S,  > 


If  the  curvature  of  the  surface  Is  uniform  then. 


d“(P,) 

n 1 n ^ 


x3/2  2-^/^  6-^/^ab/k)^/^l  - sln^ 


, ,2.  LK.  LK  .2  . LK. 
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If  as  in  the  present  Instance,  however,  the  surface  curvature  Is  nonuniform 
then. 


,3/2  3-1/2  , ,2  ^,1/4 

(21) 

where  9^^  roots  of  Equation  (17)  at  the  incident  and 

exit  points,  respectively,  of  the  curved  surface  in  the  shadow  zone;  , 

^1  (1  ■ 1,2)  being  the  corresponding  Airy  functions  and  their  derivatives 
evaluated  at  these  roots,  respectively. 

Invoking  reciprocity  for  an  individual  receiver  element  mounted  on  the 
flat  face  at  the  head  of  the  capped  cylinder  Insonifled  by  a point  source 
in  the  far-field,  we  have. 


S - 00  , - hjj  , 21  I cot  0l,  R-h^(l-jj^  sin  6^)  , 

(22) 

(r  , 6 ) being  the  coordinates  of  the  element  relative  to  the  axes  through 

o o 

the  center  of  the  ellipsoid.  Via  reciprocity. 


ikS 


11m 

S-H» 


I < 0 < ir 


(23) 


The  directivity  pattern  of  the  receiver  element  due  to  insonlficatlon  in 
the  shadow  zone,  as  defined  by  Equation  (11),  becomes 


D(0) 


cot  0 


1 1/2 


r 0-- 

exp{ikh.  (1  - sin  0 ) + Ika  j ^(1  - <^  sin^ij)) 
^ ® in 


a - ^ ,ln 


X P dJ;*^(P  )D^K(P  )exp{-e^''/^kb2/6a)^'M  4=^===^ 

t-O  " ^ “ 2 _^2  ^^^2^ 


(24) 
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the  product  D^(P,)  D^(P-)  being  defined  by  Equations  (20)  and  (21).  It 
should  be  noted  however,  that  when  we  are  dealing  with  an  edge  element  the 
terms  in  Equation  (24)  must  be  multiplied  by  the  coefficients  defined  in 
Appendix  C. 

In  order  test  the  analytical  consistency  of  Equation  (24)  we  now 
consider  the  case  where  the  insonifying  source  is  located  on  the  cylindrical 
surface  of  the  body  at  a distance  from  the  origin  of  the  coordinate 

system  along  the  axis.  This  geometry  requires  that 


.LK/ 


S " X , a 
o o 


^ - \ • (25) 


Hence,  for  a uniform  Impedance  surface,  substitution  of  Equations  (25), 


(20),  and  (16)  in  Equation  (11)  with  2^^^  6^^^  - 2^^^  gives. 


Trl/^ab/k)^/^  + *o  + ^ J 


Tr/2 


(1  - sin^  dij;)  - iTr/12 


4 X 3^/3  , ^5/6 


A 


b o 


N 

X Z 
n“0 


t tn/6,,^2,,  .1/3  LK 
expte  (kb  /6a)  q 

n 


LK  ,2  . LK. 

% *lk‘S.  > 


(26) 


If  |q^/6^^^|=  jt^l  * •|■[3'Ir(4n  + 3)/4]^^^  , then  Equation  (26)  is  identical 

8 


to  that  derived  by  Sachs.  However,  as  shown  in  Appendix  B,  these  roots 
are  applicable  when  the  baffle  is  acoustically  soft.  For  the  hard  baffle 
condition  therefore,  the  correct  roots  are  given  by  Equation  (B6) , i.e., 
|q^/6^^^|=  |t^1  * j[3Tr(4n  + l)/4]^^^  , the  erronious  choice  of  roots  in 

g 

Sachs'  paper  resulting,  as  mentioned  in  the  Introduction,  from  a failure 
to  distinguish  the  electromagnetic  from  the  acoustic  boundary  conditions. 

Returning  to  Equation  (24)  we  now  make  use  of  Equations  (10)  and  (21) 
to  compute  the  diffraction  pattern  of  a centrally  located  transducer  array 


element,  for  a uniform  Impedance  baffle  with  Z ■ Z 

pa 


varying  from  2 to  20  as  kp°  varies  from  2 to  8 • From  the  results 
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thus  obtained,  as  depicted  In  Figures  3a-3c,  which  Included  five  creeping 
wave  modes.  It  can  be  seen  that  an  order-o£-magnltude  reduction  In  the  surface 
Impedance  produces  no  reduction  In  the  back  response  of  the  element,  provided 
the  curved  surface  Is  still  effectively  rigid. 

Again,  for  a nonuniform  Impedance  surface  with  = 20  and 
varying  from  0.1  to  10,  we  obtained  the  results  depicted  In  Figures  4a-4c 
for  the  same  variation  In  kp^°  from  2 to  8 via  Equations  (10) , (21) , and 
(24).  Comparing  these  results  with  those  of  Figure  3,  It  should  be  noted 
that  the  diffraction  patterns  have  lower  main  lobe  to  side  lobe  ratios  than 
the  latter  as  a consequence  of  the  fact  that  they  are  all  normalized  with 
respect  to  the  boreslght  response  In  a relatively  rigid  baffle  (l.e., 

Zp  “ 20).  However,  the  dramatically  reduced  back  response  of  Figure  4c  Is 
a direct  consequence  of  the  almost  acoustically  soft  curved  surface. 

Finally,  although  It  can  be  seen  from  Equations  (16)  and  (17)  that 
the  attenuation  coefficients  remain  constant  for  fixed  values  of  kp^°  , 
Equations  (20)  - (22)  show  that  the  receiver  element  directivity  patterns 
will  change  via  the  diffraction  coefficients  and  differential  geometry  of 
the  surface  as  the  dimensions  of  the  ellipsoid  vary. 


Conclusions 

We  have  reviewed  the  basic  steps  required  to  determine  the  role  of 
diffraction  on  the  back  response  of  transducer  array  elements  In  uniform 
and  nonuniform  Impedance  baffles  via  the  Geometrical  Theory  of  Diffraction. 

In  particular,  we  have  considered  the  case  of  an  elllpsoldally  capped 
cylindrical  baffle,  neglecting  the  role  of  the  transition  zone  between  the 
Insonlfled  region  of  the  array  and  the  deep  shadow  zone,  which  has  yet  to 
be  deduced  analytically  for  an  arbitrary  Impedance  surface  - a task  J 
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Fig.  4c 

JimiBoimm  of  a centrelly  located  element  In  a nonuniform  impedance  baffle 
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considerably  beyond  the  scope  of  this  preliminary  Investigation.  We  have 

also  neglected  to  account  for  rapid  changes  in  the  radlus-of-curvature  on 

the  diffraction  and  attenuation  coefficients  via  'correction  factors'  such 

9 

as  those  derived  by  Voltmer  for  infinitely  hard  and  infinitely  soft  surfaces. 

As  In  the  case  of  the  transition  zone  function,  such  correction  factors  have 

yet  to  be  derived  analytically  for  an  arbitrary  impedance  surface.  Other 

implicit  limitations  of  the  investigation  also  mitigated  by  its  scope  include 

the  assumptions  that  the  surface  of  the  insonlfied  body  is  locally  reactive 

and  impenetrable,  l.e.,  the  effect  of  flexural  waves  and  Internal  sound  waves 

on  the  array  element  response  function  Is  neglected.  Nevertheless,  despite 

these  limitations,  we  believe  that  the  approach  taken  In  this  paper  Is 

capable  of  providing  a reasonable  approximation  of  an  array  element  response 

function  on  a nonuniform  Impedance,  varlable-radlus-of-curvature  surface 
o 

provided  _>  2 . 
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APPENDIX  A 

In  order  to  familiarize  the  reader  with  the  derivation  of  'diffraction' 
and  'attenuation'  coefficients  for  the  problem  of  the  elllptically  capped 
cylinder  considered  in  this  paper,  we  now  review  the  underlying  canonical 
problem  of  plane  wave  diffraction  by  an  infinite  circular  cylinder  of 
arbitrary  Impedance  z^  . Locating  an  x,  y coordinate  system  at  the 
center  of  the  cylinder,  a plane  wave  incident  field  propagating 
along  the  x axis  can  be  expressed  in  terms  of  the  polar  coordinates  (r  , 6) 


p - exp  (-ikx) 

■ exp(-lkr  cos  (J)) 

0» 

■ Z 1 ” J^(kr)  exp  (lout))  , (Al) 


where  the  time  dependence  implicit  in  exp (jut)  is  Included  post  factum, 

and  is  considered  to  be  normalized  with  respect  to  its  initial  value. 

Since  the  scattered  field  p is  composed  of  only  outward  going  waves  it 

can  be  expressed  as  a superposition  of  cylindrical  waves  weighted  by  the 

undetermined  coefficients  a as. 

m 


.-m  - 
1 a 


a (kr)  exp(lm(|>) 
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£ i **{j  (kr)  + a H^^^(kr)}  exp(lmi^) 

_ B B B 


In  this  Instance,  the  boundary  condition  at  the  surface  of  the  cylinder 

obtained  from  the  linear  constltultlve  eqixation  (l.e. , p • z v)  assumes 

a 

the  form. 


_a  -E. 

Ik  3r 


z /p  c 
a o o 


Hence, 


f u.(x)  ^ 1 

+ 2_hP*'(y)^ 


Substituting  Equation  (A5)  In  Equation  (A3)  and  applying  Watson's 
transformation  (Al)  the  field  p can  be  re-expressed  via  the  Identity 

as,  / 


_^pi€  r(„  k“ft)  ^ 

Uie  1 . L"  1H«)(X)  + 


l<»(X)^y^»'(X)  1 

mj^'cx) +z/i^’'(x)  _ 


((.l)-v  elv*  + (.i)V  _ 


Since  the  Hankel  function  of  the  first  kind  Is  analytic  In  the  complex 

v-plane  for  large  positive  values  of  kr  , It  contributes  nothing  to  the 

Integral.  However,  the  second  term  In  the  Integrand  has  simple  poles 

where  IH^^(X)  + Z (X)  I . ■ 0 . Denoting  the  roots  of  this 

V a V I X“ka 


equation  by  , the  residues  of  the  Integrand  can  be  evaluated  by 
expanding  the  denominator  In  a Taylor  series  about  such  that 
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- + Z 


' V ^ + z^h<2>'(x)}  + . . . 


(V  - V ) (X)  + z ' (X)  } (A7) 

X-ka 


For  ka  > 1 since  all  values  of  have  large  negative  Imaginary  components, 

1 - e ^ 1 . Making  this  approximation,  and  utilizing  the 

Wronsklan  relationship  ^ (ka)  H^^^'(ka)  - ^ ' (ka)  H^^^(ka)  - 4/lvka  , 

'^n  n n n 

Equation  (A6)  can  be  reexpressed  for  ka  > 1 as 


± : „(2).  , ,-lv„(l./2  - t)  ^ ,-lv„C./2  * ♦) 

*“  n-l  '’n  ’ (ka)  (ka)  1 - ll<^>  (U)  ' (ka)  1 


This  solution  can  be  further  simplified  by  redefining  the  Hankel  functions 
and  their  derivatives  In  terms  of  Airy  functions  as  follows: 


■|jH^^^ka)  - -H^2^’(ka)  - 2(2/ka)2/^  A^(-q^)  (A9a) 

V“Vn  n 


'--kTV-V 

n 

H^^^ka)  - 2(2/ka)^''^ 

n 


(A9b) 


(A9c) 


and  H^^^'(ka)  - -2(2/ka)^'^  A!(-<j  ) 

n “ 

where  A^(-x)  • ^ I cos(t^3  - xt)  dt  , for  x real  . 

/ft 


(A9d) 


(AlO) 
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Substituting  Equations  (A9a)  - (A9d)  In  Equation  (A8)  thus  gives 


/onI/3  1-2TT/3  - 
- - -zOraII)  e i 


n 


. {e-^''n(^/2  - 4.)  ^-lv^(7r/2+^: 


-Ikr 


-1/2  ,-5/6  1/3  -nr/12 
ir  2 a e 


1/6 


00 

£ ^ 

n«l 


^-lv^(ir/2  - (|»)  ^ ^-iv^(Tr/2  + 4>) 


n 4^ 


+ A:(-q^) 


-Ikr 


e Z D ^ (^-iVn(ir/2-(|»)  ^ ^-lv„(Tr/2  + <|))j 

/T  n-1  “ 


(Allb) 


where  the  'diffraction  coefficients'  D are  defined  as 

n 


2 ,r-l«  2-5^^ 


(A12) 


In  addition. 


V ■ ka  - la  a , 
n n 


CA13) 


where  the  'attenuation  coefficients'  are  defined  as 


a - — (ka/2)^^^  e^^^® 

n a 


(A14) 


Again,  from  Equations  (A9c)  and  (A9d),  since  IH^  ^(ka)  Z ^ (ka)  * 0 , 

n n 

It  follows  that 


(Alla] 
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L 


2 - -1 


n 


' (ka) 
n 


(A15a) 


.11^/6 


‘I'-V 


<2/l.)l»  *1<-V  ■ 


(A15b) 


Alternatively , 


^-1 

a 


(A16) 


Now  Equations  (A12)  and  (A14)  have  exactly  the  same  form  as  the  respective 


diffraction  and  attenuation  coefficients  summarized  by  Pathak  and 


Kouyoumjlan.^  On  the  other  hand.  Levy  and  Keller^  employed  a slightly 


different  definition  of  the  Airy  function  than  that  of  Equation  (AlO) , 

l.e. , 


■'o 


(A17) 


Relating  the  terms  In  Equations  (AlO)  and  (A17)  It  follows  that 


“ 3I/3  ^ 


(A18) 


_ — !I A • f_Y 

^2/3  ^1^  ^LK'^  ^ ’ 


(A19) 


Tjr  Ttr 

the  roots  q^  of  ) ■ 0 being  related  to  the  roots  of 


A.(-B_)  - 0 as 
1 n 


Si 


LK  . 3I/3 


Si  • 


(A20) 
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I 

9. 

ft 


V 

I 


2 

Levy  and  Keller's  diffraction  coefficients  can  thus  be  expressed  as 


.11^  ,3/2  ;-l/2  6-1/3  ^-1/6  ^1/3  ^-«/12 


,,,2,  _IK,  . LK  .2  , UL 

“Ut'-S.  > * %,  *lX<-'>n  > 


(A21) 


and  their  attenuation  coefficients  as 


LK 

Si 


. :^(ka/6)l/3  ^"x/6  ^ 


(A22) 


the  impedance  relationship  being 


a 


(A23) 


These  diffraction  and  attenuation  coefficients  can  then  be  used  locally 
on  a variable  radlus-of-curvature  surface. 
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APPENDIX  B 


In  order  to  obtain  the  roots  of  Equation  (A23)  we  Instead  employ  the 
more  fundamental  relationship  Equation  (A15a)  from  which  It  was  derived,  l.e.. 


H^^^'(ka) 


IZ 


-1 


ka  + (ka)^^^  T 


(Bl) 


where  by  comparison  with  Equations  (A13) , (A14) , and  (A20) , as  deduced  by 
Keller,^ 


(B2) 


Following  Keller,'*'  It  can  also  be  shown  that  for  ka  » 1 


„(2)..  . 2^^^  rTi  1,,  ,3/2, 

\ (ka)  - 1/2.  ,1/3  ^ 1/4  ‘^°®^4  " 3^^V  ^ 

n IT  Ckaj  T 

n 


(B3) 


and 


3X 


,2^/^  e^^^^  ,1/2,  . ir  1..  ,3/2, 

V-V  Tr^'2(ka)^'^  T-  n 4 3 n 

n n 

X-ka 


. 1 


Substituting  Equations  (B3)  and  (B4)  In  Equation  (Bl)  thus  gives  Keller's 
'tangent  approximation*  for  large  ka  , 


(2t^)^^^  tan{J  - f(2Tj^)^^^>  " -(ka)^^^  Z~^  (B5) 

or  (2t  tan{f(2T  J } - (ka)^/^  2“^  . 

u j n 4 a 

This  equation  can  then  be  solved  numerically  for  as  shown  by  Bremmer. 


Bl 
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For  the  uynptotlc  limit  of  an  acoustically  rigid  surface  (l.e.,  Z * <»), 

A 

Equation  (B5)  becomes 


slntf(2T^)^'2  - f]  - 0 . 


Hence,  ^ 


giving  - ^3ir(n  ^ 


.1/3 


n - 0,  1,  2, 


(B6) 


Thus,  <1^  ■ {3^(0  + from  Equation  (B2)  for  ■ » . 

Alternatively,  for  the  asymptotic  limit  of  an  acoustically  soft  surface 

(l.e.,  Z * 0),  Equation  (BS)  becomes 

a 

cos[  - f]  " 0 . 


Hence, 

giving 

Thus, 


^3/2  TT  , 1. 

3^^V  - 4 - (n  + 2)  ^ 

lf»  / ^ 3.  ,2/3  -iir/3 

+ j)}  e , 

.1/3 


n ~ 0,  1,  2,  . . 


(B7) 


■ 4 {3'n(n  + from  Equation  (B2)  for  Z ■ 


It  should  be  noted  that  the  values  of  and  appearing  in  Table  lb 

B1 

are  exact  estimates  obtained  by  Bremmer  from  the  Airy  function  Equation  (A23) . 
Those  obtained  from  Equations  (B6)  and  (B7)  are  approximate  to  a few  signifi- 
cant figures. 
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APPENDIX  C 

For  the  case  of  an  edge  element,  or  an  element  located  in  close  proximity 
to  an  edge,  the  singularity  in  Equation  (24)  which  occurs  as  the  angle  sub- 
tended by  the  element  at  the  center  of  the  coordinate  system  9^  approaches 
sin  removed  If  each  term  In  the  creeping  wave  series  Is 

multiplied  by  Levy  and  Keller's^  correction  factors  R^, 

IkR 

where  R • - 

n 

AT 

where 


with 


Kkp*,  q^)  » {e"^  /^2(kp*/2)‘^/^Ajj^(-q^)[l  - (kp*/2)"^^ 


+ e^^^^^(kp*/2)"^/^Aj^(-q^)[l  + (kp*/2)“^^V^’^^^]}  (C2) 


P*  • (a^/b)  (1  - sln^ 


(C3) 
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